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Given the approaching turn-on of the 
LHC, the issue of precision predictions for the 
effects of multiple gluon and multiple pho- 
ton radiative processes is a more immediate. 
A problem of some interest is the construc- 
tion of an event-by-event precision theory of 
potential luminosity processes such as sin- 
gle heavy gauge boson production processes. 
Presuming the LHC luminosity experimental 
error to reach 2% [1], the attendant theoret- 
ical precision tag on processes such as single 
W,Z production should be § — 1% so that the 
theoretical error does not adversely affect the 
respective luminosity determination and the 
attendant precision LHC physics a . Thus, we 
have developed [3] a theory of the simultane- 
ous resummation of multiple gluon and mul- 
tiple photon radiative effects, QED ® QCD 
exponentiation, to realize systematically the 
needed higher order corrections on an event- 
by-event basis, in the presence of parton 
showers, to the desired accuracy. 

We note that the new QED ® QCD ex- 
ponentiation theory is an exact resummation 
theory in the spirit of the original YFS ex- 
ponentiation for QED and is an extension 
of the QCD exponentiation theory presented 
in Refs. [4]. One should compare the lat- 
ter to the formal proof of exponentiation in 
non-Abelian gauge theories in the eikonal ap- 



a We note for reference that the current precision tag 
for the luminosity at FNAL is at the ~ 7% level [2] . 



proximation as given in Rcf. [5]: the results 
in Rcf. [4] are in contrast exact but have an 
exponent that only contains the leading con- 
tribution of the exponent in Ref. [5]- if de- 
sired, all of the contributions of the latter ex- 
ponent can be incorporated into the exponent 
Rcf. [4] while maintaining exactness. Our 
precision requirement for the LHC luminos- 
ity processes requires exact 0(a 2 , a s a, a 2 ) 
QED <E> QCD exponentiation in the pres- 
ence of parton showers, on an event-by-event 
basis, realized by MC methods for realistic 
multiple gluon and multiple photon radia- 
tive effects. Development of such MC's is in 
progress. 

We note that in Ref. [3] we have shown 
that for a process such as q + q' — > V + n(G) + 
n'(i) +X -> II' + n(g) + n'{j) + X we have 
the result 
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where the new YFS [6, 7] residuals, defined 
in Rcf. [3], $ n ,m(ki, ■ ■ ■ , k n ] k[, . . . , k' m ), with 
n hard gluons and m hard photons, repre- 
sent the successive application of the YFS 
expansion first for QCD and subsequently for 
QED. The functions SUM IR (QCED), £> QC ed 
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are given in Rcf. [3]. The residuals 
Pn,m(ki, ■ ■ ■ , k n ; k[, . . . , k' m ) are free of all in- 
frared singularities and the result in (1) is 
a representation that is exact and that can 
therefore be used to make contact with par- 
ton shower MC's without double counting or 
the unnecessary averaging of effects such as 
the gluon azimuthal angular distribution rel- 
ative to its parent's momentum direction. 

For our prototypical processes, we have 
the standard formula 

da exp (pp -> V + X -»• W + X') = 
J dx i dxjF i (x i )Fj(xj)dacx P (x i Xjs), (2) 

id 

in the notation of Rcf. [3]. We will illus- 
trate this formula with semi-analytical meth- 
ods and the case of single Z production with 
the structure functions {Fi} from Ref. [8] for 
definitencss. A MC realization will appear 
elsewhere [9]. 

As we have explained in Ref. [10], from 
the standpoint of using the result (1) with the 
structure functions {-Fi}, we intend two pos- 
sible approaches to the implied shower /ME 
matching [11]: one is based on px match- 
ing and one based on the shower-subtracted 

residuals f3„_ m . The respective shower can be 
that from PYTHIA [12], HERWIG [13], or 
the new shower algorithms in Ref. [14]. This 
combination of showers and exact ME's can 
be improved systematically to arbitrary pre- 
cision order- by-order in (a s , a) in the pres- 
ence of exact phase space [9] . 

To illustrate the size of the interplay be- 
tween QED and QCD in the threshold region, 
we compute, with and without the QED, the 
respective ratio r exp = <J exp / 'a Born, with the 
results 



' exp 




result 
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dv^fQCEDV 
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((1 - v)xix 2 s) (4) 

where we intend the well-known results 
for the respective parton-level Born cross 
sections and the value of v max implied 
by the experimental cuts under study. 
The value for the QED(g>QCD exponent 

is Iqced = {2Q}f +2C F ^L nls where 
L n i s = \nx\X2S / [i 2 when ^ is the factoriza- 
tion scale. The functions Fyfs{iqced) and 
Syfs{iqced) are well-known [7] as well: 



Fyfs(iqced) 



g— IQCEDlE 

r(i + iqced) ' 



(5) 



, QCED = QCD+QED, LHC 
,QCD, LHC 

, QCED = QCD+QED, Tevatron 
,QCD, Tevatron. 

(3) 

Here, we have used from (1) the /?o,o-level 



with Syfs(iqced) = jJqced t w/: 
C F ^)(2C(2) - i), where C(2) is Riemann's 
zeta function of argument 2, i.e., 7r 2 /6, and 
is Euler's constant, i.e., 0.5772... . 
We see that the size of the QED effects, 
~ 0.3%, is appreciable in a 1% precision 
tag budget and these effects are compara- 
ble to the QED effects found for the struc- 
ture function evolution itself in Refs. [15-19]. 
We also see that our soft, QCD correction 
is entirely consistent with the exact results 
in Ref. [20-22] and that our QED effects 
are consistent with the exact 0(a) results in 
Refs. [23-25]. Such cross-checks are essential 
in establishing a rigorous precision tag on the 
theoretical predictions. 

With an eye toward the 1% precision tag 
for the single heavy gauge boson production 
at the LHC, we have analyzed [26] the IR 
limit of the DGLAP [27] kernels themselves. 
The effects which we address are illustrated 
in Fig. 1. We apply the QCD exponentiation 
master formula from Ref. [4], the analog of 
(1) for just QCD, to the gluon emission tran- 
sition that corresponds to P qq (z), i.e., to the 
squared amplitude for q — > q{z) + G(l — z); 
as we show in Ref. [26], this allows us toget 
the following IR-improved DGLAP kernels 



G f F Y fs(,1' 
/g(7<i)*(l 



q ,e* s <{- 



(1 



>■)}, 



(6) 
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Figure 1. In (a), we show the usual process q — > <j(l — 
z) + G(z); in (b), we show its multiple gluon improve- 
ment q -» g(l-z)+Gi(£i)-|- • •+G n (e n ), z = 



^W 2 ) ' = c F fy FS ( 7g )e = 



^7, 



(7) 



Pgg{z) = 2C G F YF s(lG)e^ G {^z^ 

+ _±_(l_ 2 )7 G+ l (z l+7 G(1 _ z) 
+ Z (l-z) 1+ ^)-/ G ( 7G )5(l-z)}, 



P^( z )=F yFS ( 7G )e^i{z 2 (l-z) 
+ (l-z)V G }, 



(8) 



(9) 



where from the standard DGLAP method- 
ology [26] we have f q ( 7q ) = ± - ^ + 

7G = C G ^t=^,5 G = ^ + ^(^-I), 
and 



/g(tg) = 
+ 



l 



C G (1 + 7g)(2 + 7g)(3 + 7g) 
2 



7g(1 + 7g)(2 + 7g) 
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+ 



(1 + 7G )(2 + 7g) 2(3 + 7 g)(4 + 7g) 
1 



(2 + 7G )(3 + 7 G )(4 + 7G )' 



(10) 



Here, /3rj = 11 — \ n h where n/ is the 
number of active quark flavors. We show 
in Ref [26, 28] that these improved kernels 
change significantly the evolution of the mo- 
ments of the structure functions. For ex- 
ample, the non-singlet anomalous dimension 



A^ s becomes, from the exponentiated ker- 
nels, 

An S = C F F YFS { lq )e^ 

[B(n, lq )+B{n + 2, lq )-f q { lq )] (11) 

where B(x,y) is the beta function given by 
B(x,y) = r(x)T(y)/r{x + y). Compare the 
usual result 



A^°=C F [-\ 



n(n + 1) 



- 1 

2E ]• ( 12 ) 



i=2 ' 



The IR-improved n-th moment goes for 
large n to a multiple of —fq, consistent with 
lim^ooz™- 1 = for < z < 1; the 
usual result diverges as — 2C F lnn. The 
two results differ for finite n as well: we 
get, for example, for a s = .118, A% s = 
C F (-1.33), C F (-0.966) for (12) and (11), re- 
spectively. See Refs. [26,28] for further dis- 
cussion. 

Contact with the exact 2-loop and 3- 
loop results in Refs. [29, 30] has also been 
made [26]: for the non-singlet case, we have , 
in the notation of Ref. [30] , the exact 3-loop 
IR-improved result 



{{^) 2 {{l-z)^P^ + {z)+B 2 8{l-z)} 

+(g) 3 {(l - z)^Pif + (z) + B 3 S(l - *)}] 

(13) 

where P^ p (z) is given in (9) and the re- 
summed residuals Pns + , i = 1,2 are related 

1,2, 



to the exact results [29,30] for P„s + , i 



via 

6(0 + 



(z) = P,[ l ] + (z) - B 1+i S(l -z) + Alt' (z) 

(14) 

where the B, [29,30] and the A { nl + (z) arc 
given explicitly in Ref. [26]. The detailed 
phcnomcnological consequences of the fully 
exponentiated 2- and 3-loop DGLAP kernel 
set will appear elsewhere [9]. 

In summary, the methods illustrated and 
presented herein are all under investigation 
and implementation as we prepare for the era 
of 1% precision LHC theoretical predictions 
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on the necessary processes. Finally, one of us 
( B.F.L.W.) thanks Prof. W. Hollik for the 
support and kind hospitality of the MPI, Mu- 
nich, while a part of this work was completed. 
We also thank Prof. S. Jadach for useful dis- 
cussions. This work is partly supported by 
US DOE grant DE-FG02-05ER41399 and by 
NATO grant PST.CLG.980342. 
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